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Abstract 


After observing and investigating the double-exotic atom with the ground state lifetime T 

of about 3 X 10^*^ s, the upgraded DIRAC experiment at the CERN PS accelerator observes for 
the first time long-lived states of the same atom with lifetimes of about 10^" s and more. The 
number of characteristic pion pairs resulting from the breakup (ionisation) of long-lived atoms 

amounts to 436 ± 61, corresponding to a signal-to-error ratio of better than 7 standard deviations. This 
observation opens a new possibility to measure energy differences between p and s atomic states and 
so to determine nn scattering lengths. 


(To be submitted to Physics Letters B) 
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1 Introduction 

The experimental program of the DIRAC collaboration comprises the observation and detailed analysis 
of dimesonic atoms, which are produced by protons interacting with target nuclei O]. First, DIRAC has 
investigated atoms (pionium, A 2 ti) and measured their lifetime, T = (3.15 ^q 26) fs, in the ground 

state and hence a combination of nil scattering lengths O. Second, evidence for the production of 
dimesonic atoms with strangeness, i.e. TiK atoms, has been found Q. The data analysis has yielded a 
first measurement of the nK atom lifetime and tlK scattering lengths |@1. Third, a search for long-lived 
TTTT atom states has been performed. DIRAC plans to study the Lamb shift in Tin atoms and then to 
extract another nn scattering length combination. In this paper DIRAC presents the first observation of 
more than 400 long-lived n^n^ atoms. 

The decay probability of atoms is dominated by the annihilation process 

n^ + n^^n^ + n^ ( 1 ) 

(branching ratio ~ 99.6%) and depends on the difference between the S-wave nn scattering lengths with 
isospins zero (ao) and two (a 2 ) EIH |7l IH : 

^^Wj,o^o=R\ao-a 2 \^ with Roc |(0)^. (2) 

The expression (0) is the pure Coulomb atomic wave function at the origin with principal quantum 
number n and angular momentum quantum number 1. The most accurate ratio R has been derived with a 
precision of 1.2% in f9i. 

In order to get values of ao and a 2 separately from n^n^ atom data, one may exploit the fact that the 
energy splitting between the levels ns and np, A£'”®^"p = Ens —E„p, depends on another combination 
of the scattering lengths: 2ao + ^2 ifTOl - The influence of sfrong and elecfromagnefic interactions on the 
A 2 ;r energy structure has been studied in ifTOl [TTl [T^ [T3l IT4l and a detailed analysis performed in lUSll . 
The energy shift for the levels with the principal quantum number n and orbital quantum number I is 
composed of three contributions: AEni = AE®™ -|-AE™ + AE^f. The term AE^™ includes relativistic 
insertions, finite-size effect, self-energy corrections due to transverse photons and transverse photon 
exchange. Contributions from vacuum polarization are covered by AE™. The last term AE®]’’ takes 
into account effects from strong interaction and is related to the nn scattering lengths according to: 
AE*q =A„(2flo + « 2 )- The theoretical value for the 2s — 2p energy splitting amounts to AE^®^^p = 
(—0.59 ±0.01) eV, whereas the splitting AE”®~”p for higher principal quantum numbers n decreases 
m. By measuring the value of AE”® "p, one can determine the numerical value of AE®q, as all other 
terms in AE„/ have been calculated with high precision: the strong interaction effects contribute up to 
80% of the full energy shift. This fact provides a high sensitivity of a AE^'^^^p measurement to the value 
of the term 2ao ± a 2 - Thus, a measurement of the energy shift AE”®^''p allows to obtain a value for the 
new combination of scattering lengths 2ao ±<32- 

A method to measure AE”®^'’p has been proposed in |[T1. By investigating the influence of an applied 
electric or magnetic field on the decay probability of long-lived 71+atom states (f > 1), it is 
possible to extract an experimental value for AE"*^"p 

In inclusive processes A 2 ;i: atoms are produced in ^-states distributed over the principal quantum number 
n according to When moving inside the target, relativistic A 2 ;r interact with the electric field of target 
atoms, and some of them (N^) leave the target with an orbital quantum number Z > 0. The main primary 
excitation process is the transition ns — )■ n'p. Then, for excited states with Z > 1, the decays of A^^j. into 
two 71**, 71° ± 7 and two 7 are suppressed in accordance with (j^ because of | Yni (0) P = 0 HU- Therefore, 
the decay mechanism of such excited states is the radiative deexcitation to an ns state, annihilating 
subsequently with the lifetime T • n^ into two 71°. Thus, the A^^j decay probability is given by the shortest 
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radiative lifetime, the 2p lifetime T 2 p = 1-17 • 10^^^ s. For an average A 2 ;r momentum of 4.5 GeV/c 
( 7 ~ 16), the decay lengths are 5.7 cm (2p), 19 cm (3p) and 43 cm (4p). Using a ~100 /^m thick Be 
target and inserting a ~2 pm thick Pt foil downstream of this target ||2T|, a large fraction of the long- 
lived atoms generated in Be, reaches the Pt foil and breaks up, thus providing an extra number 
of atomic pairs (see Fig.[T]). 


Be target lOSfam Magnetic field Pt-foil 2.1um 



Fig. 1: Method to observe long-lived by means of a breakup foil (Pt). 

In order to be able to claim “observation of long-lived A^;^” - the goal of this investigation - the DIRAC 
setup has been modified in 2012 as described in the following section. 

The most accurate theoretical predictions for the S'-wave TiTi scattering lengths have been achieved in 
nil (in units M^}): 

ao = 0.220 ±0.005, ^2 =-0.0444 ±0.0010, ao - ^2 = 0.265 ±0.004 . (3) 

The best experimental results with a precision of around 4% have been obtained from studying the decays 
—)• ifTTl and —)• ifT^ as well as from measuring the atom lifetime f2l: 

|«o — = (0-2533lQQjQy). In the case of K decays, additional theoretical information can improve the 

experimental scattering length values ifTTlIT^ . 

2 Setup for detection of long-lived ;r+:;r atoms 

The DIRAC setup |[22l is sketched in Fig.|^ A high-resolution magnetic spectrometer {Ap/p ~ 3 • 10^^) 
is used to split oppositely charged meson pairs and and to measure their relative c.m. 

momentum (Q) with good precision in order to extract a dimeson atom signal 121141 . The 24 GeV/c 
CERN PS proton beam interacting with the target produces additionally free (unbound) “Coulomb pairs” 
from short-lived resonances, “non-Coulomb pairs” from long-lived sources and accidental coincidences 
(different proton-nucleus interactions). Therefore, the atom signal suffers from a continuum back¬ 
ground causing the main signal uncertainty. 

In the search for long-lived the primary proton beam hits a 103 pm thick 99.98% pure Be target, 
providing the needed A^^j yield at an acceptable proton beam intensity 1211 . The target radiation thickness 
amounts to 3.0 • 10^^ Xq (radiation length) and the nuclear interaction probability to = 2.5 • 10^^. 
The secondary channel with the whole setup is vertically inclined relative to the proton beam by 5.7° 
upward. By means of a rectangular beam collimator inside of the second steel shielding wall (Fig. 
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Fig. 2: General view of the DIRAC setup; 1 - target station with insertion, showing the Be target, magnetic field 
and Pt breakup foil; 2 - hrst shielding; 3 - microdrift chambers (MDC); 4 - scintillating fiber detector (SFD); 
5 - ionisation hodoscope (IH); 6 - second shielding; 7 - vacuum tube; 8 - spectrometer magnet; 9 - vacuum 
chamber; 10 - drift chambers (DC); 11 - vertical hodoscope (VH); 12 - horizontal hodoscope (HH); 13 - aerogel 
Cherenkov; 14 - heavy gas Cherenkov; 15 - nitrogen Cherenkov; 16 - preshower (PSh); 17 - muon detector. (The 
plotted symmetric and asymmetric events are a Kit and nK pair, respectively.) 


item 6), secondary particles are confined to ±1° in the horizontal (X) and vertical (Y) planes and thus 
to the solid angle = 1.2 • 10^^ sr. Downstream of the Be target, a 2.1 /im thick Pt foil (6.9 • 10^^ Xq, 
= 0.23 • 10^"^) for A 2 ;j breakup has been placed at a distance of 96 mm. The foil is installed at 7.5 mm 
above the primary proton beam to avoid interaction of the beam halo with Pt (Fig.[T]). The upper limit of 
the beam size in the vertical direction is Oy = 1.75 mm 1301 . The beam position in the vertical plane has 
been permanently monitored during the run by checking SFD and IH counting rates and by reconstructing 
the beam position with track information (Fig.[^. Between target and breakup foil, a permanent magnet 
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Fig. 3: Vertical coordinate of the beam position (dots) on the Beryllium target, reconstructed with track information 
for all data runs. The low edge of the Platinum foil is shown as dashed line. 


||23]| has been introduced to suppress significantly background 71 +tt pairs, generated in Be, in the low 
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Q region, where the atomic pairs from breakup are expected. The retractable magnet with a pole 
distance of 60 mm consists of a Samarium-Cobalt alloy (Sm 2 Coi 7 ) and has a maximum horizontal field 
strength of 0.25 T (see insertion in Fig. [^. The bending power of 0.02 Tm is in the secondary beam 
region relatively homogenous with a precision of better than ± 2%. 

In order to measure the shift of the vertical component Qy, which is enlarged by the horizontal magnetic 
field, e+e^ Dalitz pairs generafed in the Be target as well as e+e^ pairs produced in the Pt foil by 
photons have been investigated. The experimental Qy distribution (Fig. [^) shows a first peak around 
originating from pairs produced far downstream of the magnet and a second peak from Dalitz pairs 
at Qy = 12.9 MeV/c. After subtracting the left part of the central peak {Qy ~ 0) from the right part 
(mirrored subtraction), a peak at Qy = 2.3 MeV/c appears from e+e^ pairs produced in Pt and crossing 
only the fringing magnetic field f2M (Fig. |^). The experimental peak positions as well as the shapes 




Fig. 4: a) Experimental distribution of e+e pairs on Qy. b) Same distribution after mirror subtraction. 

coincide with simulation, which has been obtained by using the magnetic field map with a precision of 
better than 1% ll25]l . Each oppositely charged particle pair, generated in the Be target and hence crossing 
the magnetic field region, obtains an additional AQy of 12.9 MeV/c, whereas a pair produced in the Pt 
foil gets a smaller AQy, caused by the fringing field, of only 2.3 MeV/c. This difference in the Qy shift 
for pairs from Be and Pt allows to suppress significantly the background level for observing atomic pairs 
from long-lived 

The peak position for Dalitz pairs at Qy = 12.9 MeV/c has been used to control the magnetic field 
stability during the 6 month data taking in 2012: the field strength of the Sm 2 Coi 7 magnej^has been 
stable within a relative precision of better than 5 • 10^^ |[25l . 

With a spill duration of 450 ms the beam intensity has been (28-30)-10^° protons/spill. 

^ Under the same radiation condition the field strength of a permanent Nd-Fe-B magnet has been decreased by more than 
50% during the 2011 run 1^ . 
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3 Event reconstruction 

The event reconstruction has been performed by means of the DIRAC Tin analysis software already used 
for the analysis of the 2001-2003 data ||2l. Only events with one or two so-called DC tracks - tracks 
reconstructed only with DC hits - in each arm are processed according to the following criterions: 

1) One or two hadron tracks are identified in the DC of both arms with corresponding hits in VH, HH 
and PSh, but no signal in ChN and Mu detectors (Fig. [^. In each arm the earliest track, which initiates 
the readout procedure, is selected for further analysis. 

2) The DC tracks are extrapolated backward to the incident proton beam position on the target using the 
inverse transfer function of the DIRAC dipole magnet. This procedure provides a first approximation of 
the particle momenta and the corresponding intersection points in MDC, SFD and IH. 

3) SFD hits are searched for in the coordinate region defined by the position accuracy: a square of size 
±1 cm, that corresponds to 5a for high momentum and to 3a for low momentum particles. The two 
tracks should not use a common SFD hit in case of more than one hit in the proper region. If the number 
of hits around the two tracks is <4 in each SFD plane and <9 in all 3 SFD planes, the event is kept. These 
cuts reduce the data sample by 1/3, which is then called “low and medium background events”. A further 
adjustment is applied in order to find the best two-track combination: The momenta of the positively and 
negatively charged particle are modified to match the DC track X-coordinates and the SFD hits in the X- 
or (7-plane. In the final analysis, the combination with the best iri the SFD planes and vertex in the 
target is selected. 



0.004 0.0045 0.005 0.0055 0.006 0.0065 0.007 0.0075 0.008 
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[GeV/c^] 


Fig. 5 : Invariant K p mass distribution in the A region. 

The setup has been aligned using properties of A (A) decays |[27l . The decay particles are reconstructed 
and the corresponding invariant mass determined. The introduced permanent magnet (Fig. implies 
a bias in the T-components of the charged decay particle momenta. The bias amplitude depends on 
the decay vertex position inside the magnet. One observes a shift of the mass peak position and an in¬ 
crease in the width of the mass distribution. Therefore, data with magnetic field only allows to check the 
equality of the A and A masses. In order to check the A mass value, a special data set without magnet 
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has been collected (about 7% of the total data set). The distribution of reconstructed A masses is pre¬ 
sented in Fig. The experimental value M\ = 1.11568 ± 0.00001 GeV/c^ is in good agreement with 
the PDG value = 1.115683 ±0.000006 GeV/c^. Investigations of data collected in 2008-2010 
Il27]l have shown that the width of the simulated A distribution is narrower than the experimental one, 
due to underestimation of the uncertainty in particle momentum reconstruction. To improve simulation, 
a Gaussian smearing of reconstructed momenta has been introduced: ^^meared _ Xhis 

smearing of simulated momenta with C/ = (7 ±4) • 10^^ leads to a reconstructed A width consistent 
with experimental data (Fig. [^. The momentum resolution has been evaluated by investigating simu¬ 
lated events. For each track the reconstructed momentum between 1.5 and 8 GeV/c is compared with 
the generated one. The relative precision lies within 2.8 • 10^^ to 4.4 • 10^^. The resolutions in the 
components of the relative pair momentum Q, without taking into account target multiple scattering, are 
o{Qx) = o{Qy) = 0.44 MeV/c and o{Ql) = 0.50 MeV/c. 

4 Simulation 

A generator, called DIPGEN Il28ll . is used to simulate 71+pairs (atomic pairs. Coulomb and non- 
Coulomb pairs), which are generated in the Be target and Pt foil. For long-lived atoms the following 
quantities are evaluated: 

1. Ratio £„/m {Be) between number of long-lived atoms, which exit the Be target, and total number 

Na of produced atoms as a function of atom quantum numbers n, l,m: l,m) = enim{Be) x Na\ 

2. Atom decay probabilities in the space between Be and Pt as a function of n, I and atom momentum 

1123; 

3. Ratio e„im{Pt) between number N^'^‘ of long-lived atoms at the Pt foil entry and total number of 

produced atoms: {n,l,m) = e„i,„{Pt) xNa', 

4. Breakup probability Wbr in Pt as a function of n, I and m. 

In Table [T] the relative population of long-lived atom states as a function of n is given at the exit of the 
Be target and at the entry of the Pt foil taking into account their decay in-between. 


Table 1: Relative populations e„(Be) and e„(Pf) as a function of n (summed over 1 and m). 


n 

2 

3 

4 

5 

>2 

e„{Be) X 10^ 

2.48 ±0(10-^) 

1.54±0.01 

0.86 ±0.03 

0.56 ±0.06 

6.8±0.6 

en{Pt) X 10^ 

0.52 ±0(10-^) 

i.io±o(io^3; 

0.78 ±0.03 

0.54 ±0.06 

4.3±0.6 


The calculated number of long-lived atoms at the exit of the Be target is = (6.8 ± 0.6) • 10^^ x Na- 
After passing the gap of 96 mm, (4.3 ±0.6)% of the produced atoms enter the Pt foil. Their breakup 
probability in Pt varies in a wide range depending on the atom \nlm) state and momentum at the Pt foil 
entry: e.g. for an input 2p state with 4.5 GeV/c, it is about 0.7, while for a 5p state the breakup probability 
exceeds 0.95. As a result, the number of atomic pairs generated in Pt is = (4.0 ± 0.7) • 10^^ x Na, 
what corresponds to an averaged breakup probability of Wbr = 0.92. The magnetic field between Be and 
Pt has not been considered. Its influence is a decrease of the number of atoms at the Pt foil entry |[Tl[T9l 
and will be studied in future. The atomic (Pt), Coulomb (Be) and non-Coulomb (Be) pair distributions of 
relative momentum projections qx, qy, qt (initial, i.e. not smeared) at the breakup point as well as after 
multiple scattering have been simulated. In the next step, all pairs generated by DIPGEN are transferred 
to the GEANT-DIRAC (setup simulator) and ARIANE (reconstruction tool) programs. The distributions 
of the reconstructed values Ql and Qj of pairs from breakup in Pt as well as of Coulomb and non- 
Coulomb pairs, generated in Be, are obtained. The majority of the atomic pairs from the breakup in Pt 
has a Qt and Ql of less than 1.5 MeV/c (Pig. [^and|^. Moving in the fringing magnetic field, each 
atomic pair receives an additional AQy = 2.3 MeV/c. To get the initial value of Qy for atomic pairs from 
Pt, the reconstructed value must be reduced by this amount, leading to the initial value of the transverse 

momentum Qj = y^2j±72y^-±3^feV7c)2. 
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5 Data analysis and results 

The experimental distributions of pairs as a function of relative momentum Q components have 
been fitted with simulated distributions of atomic pairs (n^) from Pt, Coulomb {Nc) and non-Coulomb 
pairs {Nnc) from Be. The three corresponding numbers are free parameters in the fit. 

In the 2-dimensional {\Ql\-iQt) analysis, the experimental data have been analysed using simulated 2- 
dimensional distributions. For \Qi\ < 15 MeV/c and Qj <2 MeV/c, the \Qi\ projection of the exper¬ 
imental 2-dimensional distribution as well as of the three types of simulated pairs are shown in 

Fig. 1^. One observes an excess of events - above the sum of Coulomb and non-Coulomb pairs - in the 
low Qi^ region, where atomic pairs are expected. After subtracting background, there is a statistically 
significant signal of = 436 ± 57 (Fig. [^). The signal shape is described by the simulated distribu¬ 
tion of atomic pairs resulting from the long-lived atom breakup in the Pt foil. The atomic pair selection 
efficiency for differenf cufs on Qj is known from simulafion. Using fhis efficiency fhe fofal number 
of afomic pairs generafed in Pf is = 488 ± 64. Fig. presenfs fhe Qj projecfion of fhe same 
2-dimensional disfribufions for \Qi\ <2 MeV/c and Qj MeV/c wifh fhe same free fil paramefers. 
After background subfracfion in fhe low Q region, one observes again a sfafisfically significanl signal 
wifh a shape described by fhe simulated Qj disfribufion of afomic pairs from long-lived afom breakup in 
Pf (Fig.[^). The number of atomic pairs in fhe region Qr < 4 MeV/c is = 429 ±56. 

In fhe 1-dimensional analysis, fhe \Qi\ disfribufion is fiffed using differenf cufs Qj < 0.5, 1.0, 1.5, 
2.0 MeV/c fo sfudy fhe sfabilify of fhe afomic pair number for differenf background levels. The defected 
numbers of afomic pairs and fhe corresponding fofal numbers are shown in Table The background 
for Qt <2 MeV/c is 17 limes higher lhan for Qj < 0.5 MeV/c. Neverlheless, fhe values in fhe 1- and 
2-dimensional analyses coincide wifhin sfafisfics. This confirms fhe signal sfabilify for differenf Qj cufs, 
i.e. for differenf background levels. 

Table 2: Analysis of data collected in 2012 for different fitting procedures. The detected numbers of atomic 
pairs and the corresponding total numbers (via selection efficiency) are presented together with the back¬ 
ground contribution (Coulomb, non-Coulomb and accidental pairs) and the bt quality X^/n (n = degrees of free¬ 
dom). Errors are only statistical. 


Qt cut 
(MeV/c) 


L.tot 

Back¬ 

ground 


2-dimensional fit over Qi, Qj 

2.0 

436 ±57 

488 ±64 

16790 

138/140 

1 -dimensional fit over Ql 

0.5 

152 ±29 

467 ±88 

971 

29/27 

1.0 

349 ±53 

489 ±75 

3692 

19/27 

1.5 

386 ±78 

454±91 

9302 

22/27 

2.0 

442 ±105 

495 ±117 

16774 

22/27 


The measurement of the atomic pair number is affected by the description accuracy of the simulated 
distributions for Coulomb, non-Coulomb and atomic pairs. If the shapes of the simulated distributions 
differ from the experimental ones, then the fit parameter values might be biased. In the long-lived atom 
search, the shapes of Coulomb and non-Coulomb pairs from Be are similar in phase-space (Fig. and 
1^. Therefore, a systematic error might only arise from an incorrect description of the atomic pair distri¬ 
bution, due to the uncertainty in the Qi and Qj resolution 1^ . 
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Fig. 6: \Ql \ distribution of pairs for Qt < 2.0MeV/c. The plot a) shows the experimental distribution (points 
with statistical error) and the simulated background (solid line). The plot b) shows the experimental distribution 
after background subtraction (points with statistical error) and the simulated disttibution of atomic pairs (dotted- 
dashed line). The fit procedure has been applied to the 2-dimensional {\Ql\tQt) distribution. 
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Fig. 7: Qt distribution of pairs for \Qi\ < 2 MeV/c. The plot a) shows the experimental distribution (points 
with statistical error) and the simulated background (solid line). The plot b) shows the experimental distribution 
after background subtraction (points with statistical error) and the simulated distribution of atomic pairs (dotted- 
dashed line). The fit procedure has been applied to the 2-dimensional {\Ql\,Qt) distribution. 
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There are two main sources of systematic errors in the number of atomic pairs: 1) The A width 
correction accuracy (Section|^ leads to a systematic error of =4.4. 2) The accuracy of the measured 
Pt foil thickness dominates the uncertainty in Qt resolution, causing a systematic error of = 22 in 
the 2-dimensional analysis. In the 1-dimensional \ Ql\ analysis, this error is nearly 0. 

Another problem might arise from a hypothetical admixture of Coulomb pairs generated from beam halo 
protons interacting with the Pt foil. A peak induced by Coulomb final state interaction would be at the 
same place as atomic pairs from long-lived atoms. Beam halo level and interaction rate with Pt have 
been investigated ll^ with the result that the particle flux on Pt is practically negligible under working 
condition. To minimize any wrong interpretation, data have also been analysed under the assumption 
that they originate from Coulomb pairs generated in Pt and not from long-lived atoms: this hypothesis is 
statistically unlikely 1301. 

In summary, the 2-dimensional analysis results in the following number of atomic pairs from breakup 
in the Pt foil: = 436 ± 61. This corresponds to 7.1 standard deviations, taking into account statistical 
as well as systematic errors. In order to get an estimate for the expected number of atomic pairs, pion 
pairs generated in Be have also been analysed. To evaluate the initial value of Qy {Qt), the momentum 
shift AQy, due to the magnetic field, of 12.9 MeV/c has been subtracted from the reconstructed Qy- The 
2-dimensional experimental {\Ql\,Qt) distribution has been fitted by simulated distributions of atomic. 
Coulomb and non-Coulomb pairs from Be. Their corresponding numbers ua, Nc and Nnc are free fit 
parameters. The total number of produced atoms, Na = 17043 ±410, has been obtained by using 
the precise (1%) ratio between Na and Nc{Q < Qcut), the number of detected Coulomb pairs with small 
Q ||33l. Knowning from simulation (Section|^ that (4.0 ±0.7)% of the produced atoms {Na) break up 
in the Pt foil and the {\Ql\,Qt) ht selection efficiency is 0.89, one estimates a generation of 607 ±110 
atomic pairs, and this does not contradict the measured number of 436 ±61. 

6 Conclusion 

Long-lived 71+71^ atoms have been observed for the first time in a dedicated experiment performed by 
means of the adapted DIRAC setup. Double-exotic atoms are produced in ns states by 24 GeV/c 

CERN PS protons hitting a 103 /xm thick Be target. The pair analysis yields about 17000 
atoms, based on the measured number of small Q Coulomb pairs. These atoms are moving in the target 
and interacting electromagnetically with Be atoms. About 7% of them leave the target in excited long- 
lived states. At a distance of 96 mm downstream of the Be target, a 2.1 /rm thick breakup Pt foil has been 
installed. While passing through the gap between Be and Pt, some of the bound states (Table [T]), mainly 
shorter lived states, are decaying, whereas the rest enters the Pt foil and about 90% break up, generating 
71+71^ atomic pairs: 

n\ =436±57|stat±23|syst = 436±61|tot- 

This result corresponds to a 7.1a effect and does not contradict the estimated value of 607 ± 110. The 
observation of long-lived tt+tt^ bound states gives the possibility to study the Lamb shift and herewith a 
new Tin scattering length combination. 
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